Prospects for the measurement of oscillations with the 
ATLAS detector at LHC 
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Abstract. The capabilities of the ATLAS detector to measure the oscillations in proton-proton in- 
teractions at the Large Hadron Collider were evaluated. B° candidates in the D^n + and DJa^ decay 
modes from semileptonic exclusive events were fully simulated and reconstructed, using a detailed detector 
description. The sensitivity and the expected accuracy for the measurement of the oscillation frequency 
were derived from unbinned maximum likelihood amplitude fits as functions of the integrated luminosity. 
A detailed treatment of the systematic uncertainties was performed. The dependence of the measurement 
sensitivity on various parameters was also evaluated. 



1 Introduction 

The observed B® and states are linear combinations of 
two mass eigenstates, denoted here as H and L. Due to 
the non-conservation of flavour in charged weak-current 
interactions, transitions between B® and B® states occur 
with a frequency proportional to Am s = m B — ni^. 

Experimentally, these B®- B® oscillations have not yet 
been observed directly. In the Standard Model, their fre- 
quency is predicted in Ref. [1] to be between 12.0 ps _1 
and 17.6 ps _1 with 68% CL, and lower than 20 ps _1 at 
95% CL, significantly larger than the corresponding value 
Arrid in the B d - B® system. From measurements done by 
the ALEPH, DELPHI and OPAL experiments at LEP, by 
SLD at SLC, and by CDF at the Tevatron, a combined 
lower bound of Am s > 18.0 ps _1 has been established 
at 95% CL [2]. In the B° d - B° d system, the oscillations 
have been directly observed and a rather precise value 
of Arrid — 0.472 ± 0.017 ps _1 [3] has been measured. 

The values for Arrid and Am s predicted in the Stan- 
dard Model by computing the corresponding box diagrams, 
with the top-quark contribution assumed to be dominant, 
are proportional to \V t d\ 2 and \V ts \ 2 respectively. The di- 
rect determination of Vtd and Vt s from Auid and Am s 
is, however, hampered by hadronic uncertainties. These 
uncertainties partially cancel in the ratio: 
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where M B are the S-meson masses, B B are the bag pa- 
rameters, and f B are the i?-meson form factors. Using the 
experimentally-measured masses and a value for the ratio 



which can then be converted into a constraint of \Vtd\, the 
worst-measured side of the unitarity triangle. 

In this note, an evaluation of the capability of the AT- 
LAS detector to measure the B° s oscillations in proton- 
proton interactions at the Large Hadron Collider (LHC) 
is presented. Some quantities involved in the measure- 
ment are still uncertain (cross sections, shape of the back- 
ground, the final characteristics of the detector), therefore 
the dependence of the measurement sensitivity and of the 
expected accuracy on these parameters is also evaluated. 



2 Event selection 

The signal channels considered in this analysis for the 
measurement of B®- B® oscillations are B® — ► D s ir and 
B° -> £> s ai, with D s -> (f>ir followed by 4> -> K + K~. 

The event samples were generated using PYTHIA 5.7 
[4], passed then through the ATLAS full GEANT-based 
simulation program DICE (Inner Detector (ID) only) and 
reconstructed using an algorithm based on the Kalman fil- 
ter implemented in the xKalman package from the ATLAS 
reconstruction program ATRECON. The physics model 
used for simulation, the description of the detector and 
of the reconstruction program are presented in detail in 
Ref. [5]. 

In the simulation, 6-quark pairs were produced in pp- 
collisions at y/s = 14 TcV by including direct production, 
gluon splitting and flavour excitation processes for bb pro- 
duction. The fe-quark was forced to decay semileptonically 
giving a muon with transverse momentum 1 px > 6 GcV 



£ = y B B o f B o I y B B o f B o which can be computed in lat- 
tice QCD, a better constraint for Vt s /Vtd can be obtained, 



1 The coordinate system has the z direction along the beam 
axis, with x-axis pointing to the centre of the accelerator ring 
and j/-axis pointing upwards. The transverse momenta are com- 
puted with respect to the z axis. 
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and pseudo-rapidity < 2.4 which is used by the level-1 
trigger to select the B hadronic channels in ATLAS. The 
associated b was forced to produce the required i?-decay 
channels. All the charged final-state particles from the B 
decay were required to have pr > 0.5 GeV and \r]\ < 2.5. 

The reconstruction of the B® vertex proceeded via the 
following steps (charge-conjugate states are implicitly in- 
cluded). The <f> decay vertex was first reconstructed by 
considering all combinations of pairs of oppositely-charged 
tracks with pt > 1-5 GeV for both tracks. Kinematic cuts 
on the angles between the two tracks AipKK < 10° and 
AOkk < 10° were also imposed. Here (p denotes the az- 
imuthal angle and 8 the polar angle in the coordinate sys- 
tem defined previously. The two-track vertex was then fit- 
ted assigning the kaon mass to both tracks. Combinations 
passing a fit-probability [6] cut of 1% with the invariant 
mass within 3<70 of the nominal <j) mass were selected as cj> 
candidates. To all accepted 4> candidates, a third negative 
track with > 1-5 GeV from the remaining ones was 
added. The pion mass was assigned to the third track and 
a three-track vertex was refitted. Combinations of three 
tracks which had a fit probability greater than 1% and an 
invariant mass within 3er d s of the nominal D s mass were 
selected as D~ candidates. 

For each reconstructed Dj meson, a search was made 
for a-i candidates in three-particle combinations of the re- 
maining charged tracks. In a first step, p° mesons were 
reconstructed from all combinations of two tracks with 
opposite charges and with pt > 0.5 GeV for both tracks, 
each particle in the combination being assumed to be a 
pion. Kinematic cuts AO^ < 15° and Aip^ < 35° were 
used to reduce the combinatorial background. The two 
selected tracks were then fitted as originating from the 
same vertex; from the combinations passing a fit prob- 
ability cut of 1%, those with an invariant mass within 
1.5 Fbw of the nominal p° mass were selected as p° can- 
didates. Next, a positive track with pt > 0.5 GeV from 
the remaining charged tracks was added to the p° can- 
didate, assuming the pion hypothesis for the extra track. 
The three tracks were then fitted as originating from a 
common vertex, without any mass constraints. Combina- 
tions with a fit probability greater than 1% and with an 
invariant mass within 300 McV of the nominal a\ mass 
were selected as candidates. 

For the B a s — ► D~ir + channel, the B® decay vertex 
was reconstructed by considering all Dj candidates and 
adding a fourth track from the remaining tracks in the 
event. This track was required to have opposite charge 
with respect to the pion track from the D~ and pt > 
1 GeV. The four-track decay vertex was refitted includ- 
ing (f> and Dj mass constraints, and requiring that the 
total momentum of the £?|? vertex pointed to the primary 
vertex (within the primary vertex spatial resolutions of 
&x = &y = 28 /im and a z = 46 fim) and the momentum 
of D~ vertex pointed to the B° s vertex. 

For the B® — > Dja^ channel, the B® candidates were 
reconstructed combining the D~ candidates with the 
candidates. A six-track vertex fit was then performed with 
mass constraints for the tracks from <f> and D s ; due to the 



large a\ natural width, the three tracks from the a{ were 
not constrained to a\ mass. As in the B° s — > D~tt + chan- 
nel, the total momentum of the B® vertex was required 
to point to the primary vertex and the momentum of Dj 
vertex was required to point to the B° vertex. 

In order to be selected as B® candidates, the four- 
track and six-track combinations were required to give a 
probability greater than 1% for the vertex fit. The signed 
separation between the reconstructed B® vertex and the 
primary vertex, and between the Dj and B® vertex were 
required to be positive (the momentum should not point 
backward to the parent vertex). To improve the purity 
of the sample, further cuts were imposed: the accepted 
Bg candidates were required to have a proper decay time 
greater than 0.4 ps, an impact parameter smaller than 
55 /im and pt > 10 GeV. 

Background to the channels being considered for the 
measurement of Am s can come from two sources: from 
other four- or six-body £?-hadron decay channels, and from 
combinatorial background (random combinations with 
some or all particles not originating from a B decay). 
For B® — ► D~tt + , the following four-body decay chan- 
nels were considered as potential sources of background: 
B° -> D~ir+, B° -> D-n+ (with D~, Dj -> c/)tt' and 
— ► K + K~) and A° — > A+tt~ followed by A+ — > pK^n + . 
The similar six-body decay channels considered as poten- 
tial sources of background for B® — > D~af were: B^ — ► 
D-a+, B° d -> D-a+ (with D~, D~ -> ^ and <j> -> 
K + K~) and A° — > A+ir~ followed by A+ — > pK~ tt + it + k~ . 
The simulated four- and six-body background events were 
passed through the detailed detector-simulation program, 
reconstructed and analysed using the same programs, the 
same conditions and the same cuts as the signal events. 

In order to study the combinatorial background, a very 
large sample of simulated inclusive-muon events is needed. 
The results presented here are based on a sample of 1.1 
million bb —> p,X events, with pt > 6 GeV and \rj\ < 2.4 
for the muon corresponding to the trigger conditions. 

The bb — > p,X sample was analysed in the framework of 
a fast-simulation program ATLFAST++ [5], applying the 
same algorithms and the same cuts that were used for the 
fully-simulated samples. A careful check was made of the 
performance of the fast-simulation program by running it 
on signal and background samples, and comparing the re- 
sults with those from the detailed simulation. Reasonable 
agreement was obtained for the number of reconstructed 
events and the widths of the mass peaks for the recon- 
structed particles. 

A multi-level trigger is used to select the events for this 
analysis. The level-1 trigger is the inclusive muon trigger 
mentioned before. The level-2 trigger [7] reconfirms the 
muon from lcvcl-1 trigger, then in an un-guided search 
for tracks in the Inner Detector reconstructs a <j> meson 
and, adding a new track, a D s meson. The level-2 trigger 
uses dedicated online software. 63% of the signal events 
selected offline pass the level-2 trigger cuts; from the bb — > 
pX sample, (3.4 ± 0.2)% of the events are selected. The 
levcl-3 trigger (the event filter) uses a set of loose offline 
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Fig. 1. Reconstructed B a invariant-mass distribution for 
B® — > D^tv + decays. The open histogram shows the signal, 
the hatched histogram the background from B d — > D~iv + de- 
cays, and the dark histogram the fake reconstructed decays 
from the signal sample. The combinatorial background is not 
shown here. 



Fig. 2. Reconstructed B B invariant-mass distribution for 
B° — > -DjaJ decays. The open histogram shows the signal, 
the hatched histogram the background from B d — > D~a^ de- 
cays, and the dark histogram the fake reconstructed decays 
from the signal sample. The combinatorial background is not 
shown here. 



Process Cross-section Effective Simulated Rec. Rec. events 

(pb) cross-section (pb) events events for 10 fb _1 



B" -> 


12.81 


5.24 


23893 


2088 


2370 


B° d D-ir+ 


4.52 


2.08 


9989 


369 


400 


B ° d _» D-7V+ 


8.21 


3.74 


9989 


2 


3 


A^A+( p K-it+)it- 


89.44 


2.27 


9989 


1 


2 


B° -> D~a+ 


12.81 


2.08 


18784 


1512 


870 


B° d -> D~a+ 


20.98 


3.35 


9699 


190 


340 


B° d -> D-a+ 


8.20 


1.40 


9949 


1 


1 


Al -» A+{ P K~TY+TT + n-)TY- 


1.97 


0.20 


10994 








Comb. background sum 






1.1 x 10 6 


see text 


3750 



Table 1. Signal and background samples analysed for the study of B°- B° oscillations. The combinatorial background is the 
sum for both B® — > DJtt + and B® — > D^a^ analysis channels. 



cuts, reducing the rate to (0.26 -j- 0.41)% of the /iX rate, 
depending on the actual values which are set for the cuts. 

The reconstructed B® invariant-mass distributions in 
the decay channels B® — > D~n + and B® — > D~a^ are 
shown in Figure 1 and Figure 2, respectively, for an inte- 
grated luminosity of 10 fb^ 1 . 

The numbers of events expected for the various sig- 
nal and background channels that have been analysed are 
given in Table 1 for an integrated luminosity of 10 fb -1 . 
The first cross-section column in Table 1 gives the chan- 
nel cross-section, without any cuts on final-state parti- 
cles, assuming a cross-section of 2.3 fjb for the process 
bb -> n(pr > 6 GeV, |r/| < 2A)X. The effective cross- 
section is the cross-section after the cuts on charged final- 
state particles were applied during simulation. 

The events reconstructed from the samples for the 
exclusive decay modes were counted in a ±2er window 



around the nominal B® mass. Using the fraction of events 
reconstructed in the simulated sample and the number of 
events expected for an integrated luminosity of 10 fb _1 , 
the expected number of reconstructed events was esti- 
mated. Corrections for muon efficiency (on average 0.82) 
and for level-2 trigger efficiency (0.63) were also applied. 

A total of 3240 reconstructed events is expected for 
the Bg — ► D s tt and B® — * D s a\ decay channels for an 
integrated luminosity of 10 fb _ . 

The only significant background comes from the B® — > 
D~n + and B® — > D~ channels, and from the com- 
binatorial background. Note that the number of recon- 
structed events from the two decay channels is conser- 
vative since the branching-ratio values used are upper 
limits. As expected, due to the combination of the D~ 
mass shift (M D - - M Ds = 90 McV) and B° d mass shift 
(M B o - M B o w 100 MeV), very few B a d -> L>~a+, D--K+ 
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events are reconstructed in a ±2<7 window around B® 
nominal mass. Due to the different decay topology, the 
A® — > A+ir~ channel does not give any contribution to 
the background. 

The statistics available for estimating the combinato- 
rial background are very limited, despite the large size (1.1 
million events) of the fj,X sample. Each simulated event 
was therefore passed 50 times through the fast-simulation 
program, different random smearing of the track parame- 
ters being applied each time. The number of background 
events was counted in a 2a mass window around the B® 
nominal mass. On average, 0.12 events per pass were re- 
constructed in the mass window, summing the B® — > D~ n + 
and B® — ► Dja^ channels. Normalising to the number 
of /iX events expected for an integrated luminosity of 
10 fb" 1 , applying correction factors for the reconstruction 
and trigger efficiencies, the combinatorial background was 
estimated to be 3752 events, and the range of variation was 
estimated to be between 2079 and 5425 events at 90% CL. 
Correlations between the results from the 50 passes were 
taken into account. 
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Fig. 3. Decay-radius resolution for the decay channel 
B°^D-(<tnr-)-K+. 



3 Proper-time reconstruction and resolution 

The proper time of the reconstructed B® candidates was 
computed from the reconstructed transverse decay length, 
transverse momentum, p^: 

d xy M B o 



d xy , and from the B® 



t = 



cpt 



— d x y 



where g = Mgo / (cpx)- 

The transverse decay length is the distance between 
the interaction point and the &-hadron decay vertex, pro- 
jected onto the transverse plane. Figure 3 shows, for the 
example of the B° s — ► D~ (<f)Tr~)iT + decay mode, the differ- 
ence d xy — d^. y fitted with two Gaussian functions, where 
d° y is the true transverse decay length. For each event, 
the decay-length uncertainty, a d , was estimated from 
the covariance matrices of the tracks associated with the 
vertices. The pull of the transverse decay length, (d xy — 
^ X y) I a d xy ■ was found to have a Gaussian shape with a 
width of S dxy = 0.959 ± 0.017 for the B Q S -> DJ-k + chan- 
nel and S dxy = 0.954±0.020 for the B° -> D~af channel. 

The distributions for (g — go) /go also have a Gaus- 
sian shape for both B° s decay channels, with a width of 

DZtt + channel and 



(0.715 ± 0.014)% for the B° s 
(0.636 ±0.013)% for the B° 



D 



channel. Here 



go = to/d xy , with to being the true proper time. The 
(g ~ .9o ) /go distribution for the B® — ► D~ (4>Tr~)it + decay 
mode is shown in Figure 4. 

The proper-time resolution function Res(i | to) was pa- 
rametrised with the sum of two Gaussian functions, with 
parameters given in Table 2: 



Hes(t|i ) =/i 
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Fig. 4. Fractional resolution on g-factor for the decay 
channel B® — > D~ (4>tv~)tt + . 



Figures 5 and 6 show, for the decay channels B® — ► 
D~n + and B® — > D~a± , the proper-time resolution to- 
gether with the parametrization obtained from the func- 
tion given above. The distributions show deviations from 
the Gaussian shape, illustrated by the significant fraction 
of the second, broader Gaussian function in the param- 
etrization function. For B® mesons, the ratio of the two 
Gaussian functions was fixed to the ratio from B® param- 
etrization. 

The parametrization chosen above reproduces well the 
tails seen in the distribution for reconstructed events. It 
has the advantage that the integrals in the log-likelihood 
function given in Eq. (3) can be computed analytically, 
much faster than the numerical computation. The depen- 
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D7 7V+ 




D~ a+ 


Q 


Fraction f a 


Width a a 


Fraction f a 


Width a a 




(%) 


(fs) 


(%) 


(fs) 


1 


59.6 ±6.6 


51.5 ±4.0 


62.5 ± 14.1 


51.6 ± 6.4 


2 


40.4 ± 6.6 


107.3 ±8.5 


37.5 ± 14.1 


92.8 ± 12.7 




Bd — * 


D~n + 






1 


59.6 ± 1.8 


56.9 ±5.5 


62.5 ± 2.6 


60.3 ± 7.2 


2 


40.4 ± 1.8 


102.9 ±9.6 


37.5 ± 2.6 


100.1 ±31.6 



Table 2. Proper-time resolution function Res(t | to) parame- 
trization with the sum of two Gaussian functions. 



dence of the proper-time resolution on the proper time 
a(t Q ) = ^J{gS d!Dy a d;Dy ) 2 + (t S g ) 2 can be neglected in the 
first approximation, as the factor S g is very small. 
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4 Likelihood function 

The probability density to observe an initial B® meson 

(j = d, s) decaying at time to after its creation as a B Q - 
meson is given by: 



Pj{t , Mo) 



r 2 - 



2 



2 Pi 



cosh 



ArA 



where Af^ 



1 H 



-Tjto 



[Iq cos Arrijto 



r[)/2 and /x 



(1) 



-1. 



For the unmixed case (an initial B® meson decays as a B® 
meson at time to), the probability density is given by the 
above expression with yUo = ±1- The small effects of CP 
violation are neglected in the above relation. Unlike Ar,i, 
which can be safely neglected, the width difference in the 
B°- B° system Ar s could be as much as 20% of the total 
width [8]. 

The above probability is modified by experimental ef- 
fects. The probability as a function of ^o an d the recon- 
structed proper time t is obtained as the convolution of 
Pj{to, /xo) with the proper time resolution Res.,(t | to): 

Pj(t ,Ho) Res^t | t ) dt 



Fig. 5. Proper-time resolution for the decay chan- 
nel B° -> D~tt + . The curve displays the resolution 
as obtained from the Res(t | to) function given in the 
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Fig. 6. Proper-time resolution for the decay chan- 
nel B° -> D~af. The curve displays the resolution 
as obtained from the Res(t | to) function given in the 
text. 



with N a normalization factor. Assuming a fraction uij 
of wrong tags at production or decay, the probability be- 
comes: 



combinatorial background cb) in the total sample of type 
k, one obtains the probability density: 



Qj(t, M) = (1 - Wj)Qj(t, M) + Ujqj{t, -m) 



(2) 



For each signal channel, the background is composed of 
oscillating B® mesons with probability given by the ex- 
pression (2) and of combinatorial background, with prob- 
ability given by the reconstructed proper time distribu- 



E 

— s.d.c 



fj [(! ~ u^qjit, m) ± Ujqj(t, -fJ,)] 



The index k = 1 denotes the B2 



tion. For a fraction / fe of the j component (j = s, d, and k = 2 the B° — > D s channel. 



— ► D s 7r + channel and 
The likelihood of the 
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total sample is written as 

C(Am s ,Ar s ) = f[ Y[pd£ k (ti,IH) (3) 

k=l i=l 

where N^ v is the total number of events of type k, and 
iV c h = 2. Each pdf fe was properly normalized to unity. 

5 Extraction of Am s measurement limits and 
accuracy 

The maximum value of Am s measurable in ATLAS was 
estimated using a simplified Monte Carlo model. The in- 
put parameters of this model were: for each signal chan- 
nel k the number of signal events, N* lg , the number of 
background events from B® decays, Ng , and the num- 
ber of events for the combinatorial background, N* h ; the 
characteristics of the events involved in the computation 
of the proper-time resolution (see below); the wrong-tag 
fraction. The wrong-tag fraction was assumed to be the 
same for both B° and mesons: ui s = u>d = 0.22 (see 
Ref. [9]). A Monte Carlo sample with N sig = J\A + N? g 
signal events oscillating with a given frequency Am s to- 
gether with N B o = Ng + iVL background events oscil- 

« d d 

lating with frequency Am d and N cb = N\ h + N^ h combi- 
natorial events (no oscillations), was generated in the fol- 
lowing way. For each event with an oscillating b hadron, 
the true proper time was generated according to an expo- 
nential distribution with the slope [-P(l ± AT/{2r))}- 1 
which takes into account the contribution of a non-zero 
AT. Here r = r s , Ar = Ar s for B° s events and T = r d , 
AF = for B® events. The sign in the slope is chosen 
according to the contribution of each proper time, corre- 
sponding to r(l ± Ar/(2r)), to the exponential decay 
from Eq. (1), c~ Fjt ° cosh(Z\_T ? to/2), where the cosh term 
was factorized from the expression in parentheses and in- 
cluded in the exponential part. The uncertainty on the 
measurement of the transverse decay length, <7d , and the 
true value of the g-factor, go, were generated at random 
according to the distributions obtained from the simulated 
samples, fitted with appropriate combinations of Gaussian 
and exponential functions. From the computed true decay 
length, d xy = to /go, the corresponding reconstructed de- 
cay length was generated as d xy = d xy + Sd xy o'd mv O. The 
reconstructed g-factor was generated as g = go + goS g £2' . 
Both Q and fi' are random numbers distributed accord- 
ing to the normal distribution. From the transverse de- 
cay length and g-factor, the reconstructed proper time 
was then computed as t = gd xy . The probability for the 
event to be mixed or unmixed was determined from the to 
and Am s values (Anid value if the event was a B d event) 
using the expression (1 — cos(Amjt ) / cosh(Arjto/2)) /2 
left out from Eq. (1) after the exponential part is ex- 
tracted. For a fraction of the events, selected at random, 
the state was changed between mixed and unmixed, ac- 
cording to the wrong-tag fraction, Wtag- For the combi- 
natorial background, the reconstructed proper time was 



generated assuming that it has the same distribution as 
the one for J3° mesons. Half of the combinatorial events 
were added to the mixed events and half to the unmixed 
events. 



5.1 Am s measurement limits 

The Am s measurement limits were obtained applying the 
amplitude fit method [10] to the 'data sample' generated 
as above. In this method a new parameter, the B® oscilla- 
tion amplitude A, is introduced in the likelihood function 
by replacing the term '/xp cos Am s to^ with 'ixo-4cos Am s to' 
in the B® probability density function given by Eq. (1). 
For each value of Am s , the new likelihood function is 
minimized with respect to A, keeping all other param- 
eters fixed, and a value A ± <r^ at is obtained. One ex- 
pects, within the estimated uncertainty, A = 1 for Am s 
close to its true value, and A = for Am s far from 
the true value. One defines a 5<r measurement limit as 
the value of Am s for which 1 / <r_4 = 5, and a sensitivity 
at 95% confidence limit as the value of Am s for which 
^l a A = L645. Limits are computed with the statistical 
uncertainty <r^ at , and, in some cases, with the total un- 
certainty cr^ )tal = cr^ at + a b J st . The systematic uncertainty 
a s J st is described in the next section. 



5.2 Systematic uncertainties 

The systematic uncertainties on the B® oscillation ampli- 
tude a s j[ st are computed as 

A _stat 

a T * = AA+(l-A)^ (4) 
"A 

where AA is the difference between the value of the ampli- 
tude when a single parameter is changed and the analysis 
repeated, and the value for the 'nominal' set of parame- 
ters; Aa s J at is defined in a similar way. 

The following contributions to the systematic uncer- 
tainties were considered: 

— A relative error of 5% was considered for the wrong-tag 
fraction for both B® and B^. 

— The widths of the Gaussian functions from the param- 
etrization of the proper time resolution given in Table 2 
were varied by ±1ct. 

— The fraction f B o = BR(b -> B° s ), the B° s lifetime and 
the Amd value were varied separately by the uncer- 
tainty quoted in Ref. [11]. 

— An uncertainty of 5% was assumed for the decay time 
r c b of the combinatorial background. The shape re- 
mained exponential, only the decay time was modified. 

These contributions were added in quadrature to give 
the systematic uncertainty. 

Table 3 shows the dependence of the amplitude and its 
statistical and systematic uncertainties on Am s , as well 
as the contribution of each component to the systematic 
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Am s 


ps" 1 


5 ps^ 1 


10 ps^ 1 


15 ps" 1 


20 ps^ 1 


25 ps^ 1 


30 ps" 1 


35 ps^ 1 


40 ps" 1 


A 


0.0450 


-0.1146 


0.1892 


-0.0337 


0.0420 


0.0072 


-0.2906 


0.4318 


-0.5406 


stat 

a A 


±0.0476 


±0.0707 


±0.0895 


±0.1183 


±0.1673 


±0.2456 


±0.3570 


±0.5640 


±0.9300 


syst 

a A 


+0.0967 
-0.0838 


+0.1285 
-0.1018 


+0.1297 
-0.0960 


+0.1468 
-0.1144 


+0.1796 
-0.1424 


+0.2035 
-0.1580 


+0.2980 
-0.2264 


+0.3674 
-0.2658 


+0.3636 
-0.2592 


Systematic contributions 




















- wrong-tag fraction 


+0.0387 
-0.0382 


+0.0361 
-0.0360 


+0.0445 
-0.0349 


+0.0389 
-0.0391 


+0.0389 
-0.0312 


+0.0394 
-0.0329 


+0.0408 
-0.0342 


+0.0355 
-0.0319 


+0.0418 
-0.0364 


- proper time resolution 


i n nnnn 
-|-U.UUUU 

-0.0000 


-0.0127 


-hU.UOOl 

-0.0450 


-t-u.uyou 
-0.0758 


-0.1143 


-|-U.l / uu 

-0.1362 


-0.2098 


-0.2490 


-|-U.o40o 
-0.2425 


- f a fraction 


+0.0884 
-0.0742 


+0.1227 
-0.0943 


+0.1086 
-0.0773 


+0.1049 
-0.0761 


+0.1107 
-0.0790 


+0.1035 
-0.0724 


+0.1069 
-0.0780 


+0.1146 
-0.0861 


+0.0888 
-0.0661 


- B° lifetime 


+0.0038 
-0.0036 


+0.0009 
-0.0009 


+0.0034 
-0.0032 


+0.0030 
-0.0033 


+0.0060 
-0.0000 


+0.0146 
-0.0076 


+0.0010 
-0.0014 


+0.0132 
-0.0088 


+0.0372 
-0.0390 


- Am d 


+0.0055 
-0.0055 


+0.0001 
-0.0001 


+0.0003 
-0.0003 


+0.0004 
-0.0003 


+0.0002 
-0.0002 


+0.0007 
-0.0007 


+0.0001 
-0.0002 


+0.0003 
-0.0003 


+0.0020 
-0.0019 


- Tcb 


+0.0012 
-0.0013 


+0.0018 
-0.0016 


+0.0014 
-0.0021 


+0.0031 
-0.0027 


+0.0001 
-0.0002 


+0.0057 
-0.0061 


+0.0047 
-0.0041 


+0.0094 
-0.0109 


+0.0403 
-0.0339 



Table 3. The oscillation amplitude A and its statistical and systematic uncertainties as a function of Am s for an integrated 
luminosity of 10 fb" 1 . The contribution of each component to the systematic uncertainty is also given. 



uncertainty for an integrated luminosity of 10 fb~ . In 
the generated event samples, the value of Am s was set to 
Z\mf on = oo. 

From Table 3, it can be seen that the dominant con- 
tributions to the systematic uncertainty come from the 
uncertainty on f B o fraction and from the parametrization 
of the proper time resolution. 

The uncertainty on f B o fraction has the value given 
in Ref. [11]. It is, however, expected that the uncertainty 
will be much smaller at the time this analysis will be done 
with data. Even now (October 2001), the most up-to-date 
preliminary value [12] of f B o is 0.099 ± 0.011, the uncer- 
tainty being at the level of 11%, to be compared with the 
value of 17% used here for the consistency of the data. If 
one considers the intense activity in B— physics, one can 
assume that an uncertainty of ~ 5% will be achieved at 
the time of data analysis. 

The uncertainties on the parameters from the param- 
etrization of the proper time resolution depend on the 
MC statistics. A larger MC sample can be generated than 
the one used in this work, in order to reduce these un- 
certainties. It is therefore reasonable to assume that the 
uncertainty on the widths from Table 2 can be reduced 
to half of the actual values. 

The systematic uncertainties computed with these 'pro- 
jected' uncertainties are given in Table 4. For the other 
contributions, the values from Table 3 were used. 

The value of the systematic uncertainty as computed 
here should be considered with caution, as a rough esti- 
mate. 



5.3 Results 

The amplitude as a function of Am s for the nominal set 
of parameters defined in the previous sections, Ar s = 
and an integrated luminosity of 10 fb _1 is shown in Fig. 7. 



Fig. 8 shows the significance of the measurement in units 
of 04. The 5a measurement limit is 22.5 ps^ 1 and the 
95% CL sensitivity is 36.0 ps _1 , when computed with the 
statistical uncertainty only. Computed with the total un- 
certainty, the 5(7 measurement limit is 16.0 ps _1 and the 
95% CL sensitivity is 34.5 ps _1 for the actual systematic 
uncertainties, and 21 ps _1 and 35.5 ps _1 for the projected 
systematic uncertainties. 

For an integrated luminosity of 30 fb _1 , the 5a mea- 
surement limit is 29.5 ps _1 and the 95% CL sensitivity 
is 41.0 ps _1 , computed with the statistical uncertainty 
only. The 5a measurement limit and the 95% CL sensitiv- 
ity become 18.5 ps _1 and 37.5 ps _1 , respectively, for the 
actual set of systematic uncertainties and 27.0 ps _1 and 
40.5 ps _1 for the projected systematic uncertainties. 



6 Dependence of Am s measurement limits 
on experimental quantities 

Some quantities involved in the measurement of the B® 
oscillations are not yet known with enough precision. The 
cross section for bb production could be more than twice 
the assumed cross section [5]; some decay branching ra- 
tios are also not well determined. The value of Ar s is not 
yet measured. The characteristics of the detector could 
be slightly different, depending on the final configuration 
of the detector. The shape and the fraction of the back- 
ground will depend, among other factors, on the accelera- 
tor luminosity and on the characteristics of the detector. 
It is therefore necessary to estimate the dependence of the 
measurement sensitivity for various values of the param- 
eters. The limits presented in this section are computed 
with the statistical error only. 

The dependence of the Am s measurement limits on 
Ar s /r s was determined for an integrated luminosity of 
30 fb _ , other parameters having the nominal value. The 
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Am s 


ps" 1 


5 ps" 1 


10 ps" 1 


15 ps" 1 


20 ps" 1 


25 ps" 1 


30 ps" 1 


35 ps" 1 


40 ps" 1 


A 


0.0450 


-0.1146 


0.1892 


-0.0337 


0.0420 


0.0072 


-0.2906 


0.4318 


-0.5406 


stat 

°A 


±0.0476 


±0.0707 


±0.0895 


±0.1183 


±0.1673 


±0.2456 


±0.3570 


±0.5640 


±0.9300 


syst 
a A 


+0.0494 
-0.0492 


+0.0564 
-0.0471 


+0.0596 
-0.0475 


+0.0679 
-0.0617 


+0.0848 
-0.0661 


+0.0976 
-0.0823 


+0.1370 
-0.1167 


+0.1660 
-0.1411 


+0.1747 
-0.1496 


Systematic contributions 




















- proper time resolution 


+0.0000 
-0.0000 


+0.0134 
-0.0133 


+0.0229 
-0.0226 


+0.0476 
-0.0379 


+0.0681 
-0.0532 


+0.0825 
-0.0708 


+0.1277 
-0.1077 


+0.1579 
-0.1341 


+0.1586 
-0.1340 


- f a fraction 


+0.0299 
-0.0301 


+0.0412 
-0.0272 


+0.0321 
-0.0226 


+0.0285 
-0.0287 


+0.0316 
-0.0238 


+0.0305 
-0.0242 


+0.0278 
-0.0286 


+0.0333 
-0.0265 


+0.0252 
-0.0205 



Table 4. The oscillation amplitude A and its statistical and systematic uncertainties as a function of Am a for an integrated 
luminosity of 10 fb _1 , computed with the 'projected uncertainties' on proper time parametrization and on f s . See the text for 
details. 
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Fig. 7. The B° oscillation amplitude as a function of Am s for Fig. 9. The dependence of Am s measurement limits on 
an integrated luminosity of 10 fb _1 . Ar a /r s . 
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Ar s /r s was used as a fixed parameter in the amplitude fit 
method. The results are shown in Fig. 9 and the numerical 
values are given in Table 5. No sizeable effect is seen up 
to Ar s /r s < 30%, therefore AT a was set to zero for all 
other cases. 



Ar s /r s 


5<r limit 


95% CL sensitivity 


(%) 


(ps" 1 ) 


(ps" 1 ) 





29.5 


41.0 


10 


29.5 


41.0 


20 


29.5 


41.0 


30 


29.0 


41.0 


40 


29.0 


40.5 


80 


27.0 


39.0 


100 


25.5 


38.0 



so Table 5. 

Ara ' [ps "' 1 Ar s /r s . 



The dependence of Am a measurement limits on 



Fig. 8. The measurement significance as a function of Am s 
for an integrated luminosity of 10 fb . 
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Fig. 10. The dependence of Am s measurement limits on the Fig. 11. The dependence of Arris measurement limits on the 



integrated luminosity. 



The dependence of the Am s measurement limits on 
the integrated luminosity is shown in Fig. 10, with the 
numerical values given in Table 6. 



Luminosity 


5a limit 


95% CL sensitivity 


(fb- 1 ) 


(ps- 1 ) 


(ps- 1 ) 


5 


17.5 


32.0 


10 


22.5 


36.0 


20 


27.0 


39.0 


30 


29.5 


41.0 



Table 6. The dependence of Am s 
integrated luminosity. 



measurement limits on the 



The shape and the fraction of the combinatorial back- 
ground were varied within reasonable values. It was as- 
sumed that the shape remains exponential and only the 
decay time was modified. The number of events for the 
combinatorial background was varied within ±50% of the 
number determined in Section 2, keeping the number of 
B Q S and B® events to the nominal values. The dependence 
of the Am s measurement limits on the decay time of the 
combinatorial background is shown in Fig. 11, with the 
numerical values given in Table 7, while the dependence 
on the number of events is given in Fig. 12 and Table 8. 



decay time of the combinatorial background, assuming an ex- 
ponential decay. 



T cbg 


5<r limit 


95% CL sensitivity 


(PS) 


(ps" 1 ) 


(ps" 1 ) 


0.25 


31.5 


42.5 


0.50 


30.5 


42.0 


1.00 


29.5 


41.5 


1.54 


29.5 


41.0 



Table 7. The dependence of Am a measurement limits on the 
decay time of the combinatorial background, assuming an ex- 
ponential decay. 



.5 60 



<3 



40 



20 



• 5o limit 
sensitivity 



0.5 



1.5 



7 Accuracy of the Am s measurement 

For Am s values smaller than the 5ct measurement limit, 
the expected accuracy is estimated using the log-likelihood 
method, with the likelihood function given by Eq. (3). In 
the fit, the Am s value was free, while the other parame- 
ters were fixed to their nominal values. An example of the 



N cbg /Ncb g (n°m) 

Fig. 12. The dependence of Arris measurement limits on the 
number of events for the combinatorial background. 



10 B. Epp, V.M. Ghete, A. Nairz: Prospects for the measurement of B° oscillations with the ATLAS detector at LHC 





5ct limit 

(PS" 1 ) 


95% CL sensitivity 
(P^ 1 ) 


Luminosity 
(ID ) 


Amf n 
IP S ) 


Am s ±CT atat s ±cr ayst s 

\P S ) 


Obs. 


0.50 


31.0 


42.5 


5 


17.5 


17.689 ± 0.083 ± 0.002 


5f7 limit 


0.75 


30.5 


41.5 


10 


15.0 


15.021 ± 0.049 ± 0.002 




1.00 


29.5 


41.0 


22.5 


22.396 ±0.072 ±0.005 


5cr limit 


1.25 


29.0 


40.5 




15.0 


14.949 ± 0.033 ± 0.002 




1.50 


28.5 


40.0 


20 


20.0 


20.041 ±0.068 ±0.005 




Table 8. The dependence of Am s measurement limits on the 




27.0 


26.948 ± 0.070 ± 0.003 


5cr limit 


number of events for the combinatorial background. 


30 


15.0 
20.0 


14.942 ± 0.028 ± 0.004 
20.010 ± 0.043 ± 0.002 












29.5 


29.708 ± 0.083 ± 0.007 


5cr limit 




Am s [ps" 



Fig. 13. The negative log-likelihood function for an integrated 
luminosity of 10 fb _1 and a generated Z\mf on = 22.5 ps -1 . 



likelihood function is given in Figure 13 for an integrated 
luminosity of 10 fb _1 and a generated Z\mf cn = 22.5 ps _1 . 

The accuracy was determined for different values of the 
integrated luminosity. The results arc given in Table 9. 

The systematic accuracy of the Am s measurement was 
obtained similarly to the systematic uncertainty on the 
measurement limits: one parameter was changed, keeping 
the other fixed to their nominal values, and the analysis 
was repeated. The difference between the value of the fit- 
ted Am s parameter when a single parameter is changed 
and the analysis repeated, and the value for the 'nominal' 
set of parameters was taken as the systematic contribution 
corresponding to that parameter. The total systematic un- 
certainty is the quadratic sum of the individual contribu- 
tions. The parameters changed are the same as for the 
measurement limits, see Section 5.2, within the 'actual' 
uncertainties. 

One can see from Table 9 that the Am s measurement, 
if within the reach of the ATLAS detector, will be domi- 
nated by the statistical errors. 



Table 9. The accuracy of Am s measurement as a function 
of the integrated luminosity, (7 S ^ S represents the statistical 
uncertainty, a^^ s the systematic uncertainty. 



Carlo simulated events, propagated through a detailed 
simulation of the detector. For an integrated luminos- 
ity of 10 fb — , a 5<7 measurement will be possible up to 
22.5 ps _1 , the experiment being 95% CL sensitive up to 
36.0 ps _1 . These limits increase to 29.5 ps _1 and 41.0 ps _1 , 
respectively, for an integrated luminosity of 30 fb _ . The 
effect of the systematic uncertainties on the measurement 
limits is rather limited, if the precision of the f B o fraction 
is improved with respect to the actual value, and if enough 
Monte Carlo events are generated for the parametrization 
of the proper time. 

If the Am s value will be within the ATLAS reach, the 
measurement will be dominated by the statistical error. 
A total uncertainty of the order of ~ 0.07 ps _1 is ex- 
pected for a value of Am s = 22.5 ps _1 for an integrated 
luminosity of 10 fb _1 , the total uncertainty decreasing to 
~ 0.04 ps _1 for a luminosity of 30 fb _1 . 

The dependence of the measurement sensitivity on var- 
ious parameters which will be known only at the time of 
the data analysis was also evaluated. 

The values obtained in this note for the measurement 
limits and accuracies should be re-evaluated at a later 
time, taking into account the changes in the detector ge- 
ometry and in the simulation and reconstruction software. 
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